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1. 


INTRODUCTION 


This  report  describes  the  results  of  the  work  carried  out  by  Da¬ 
nish  Hydraulic  Institute  (DHI)  under  Contract  No.  N62474-84-C-3 1 4 2 
for  the  Naval  Civil  Engineering  Laboratory  (NCEL) ,  Department  of 
the  Navy,  Port  Hueneme,  U.S.A.  The  Danish  Maritime  Institute  (DMI) 
has  acted  as  consultants  to  DHI  during  the  project. 

The  scope  of  work  has  covered  the  influence  on  the  drag  force  on  a 
ship  moored  in  a  beam  current  in  shallow  water  from 

o  Blockage  effects 
o  Reynolds  number 

Furthermore,  the  possibility  of  having  instability  type  of  motions 
was  investigated.  During  the  course  of  the  project  NCEL  has  been 
kept  informed  about  the  progress  of  the  work. 


2 .  SUMMARY 

A  limited  model  test  programme  has  been  conducted  in  order  to 

establish  the  necessary  flume  width  for  measuring  lateral  drag 

forces  increased  by  less  than  5%  due  to  blocking,  as  no  blockage 

correction  formulas  for  the  three-dimensional  flow  exists.  The 

application  of  two-dimensional  correction  formulas  indicate  a  very 

conservative  estimate  of  a  flume  width  9  times  L  .  The  model 

PP 

tests  have  shown  that  for  a  flume  width,  B,  of  5  times  L  ,  late- 

P  P 

ral  drag  forces  are  not  affected  significantly  by  blockage. 

The  influence  of  Reynolds  number  on  drag  forces  is  well-known  for 
various  bluff  bodies.  The  drag  force  is  a  form  drag  with  insigni¬ 
ficant  contribution  from  shear  stresses.  For  circular  cylinders 
the  drag  coefficient  varies  significantly  with  the  Reynolds  num¬ 
ber,  with  a  subcritical  value  of  CD  equal  to  1.2  and  a  postcritial 
value  of  0.7,  approximately,  for  absolutely  smooth  cylinders.  For 
cylinders  with  a  rough  surface  the  postcritical  value  of  CD  in¬ 
creases,  approaching  the  subcritical  value.  For  bluff  bodies  with 
cross-sections  having  more  sharp  edged  corners  -  such  as  squares 
or  rectangles  -  the  drag  coefficient  varies  insignificantly  with 
the  Reynolds  number,  as  the  separation  is  fixed  to  the  sharp  corn- 


The  shape  of  a  ship  hull  with  bilge  keels  has  similarities  with 
both  of  these  types  of  bluff  bodies.  The  influence  of  the  Reynolds 
number  is  thus  a  priori  expected  to  be  much  less  than  for  a  smooth 
circular  cylinder,  and  especially  the  difference  in  the  Cp-values 
for  the  sub-  and  postcritical  Reynolds  number  is  expected  to  be 
substantially  sir.  Her.  However,  comparisons  are  extremely  diffi¬ 
cult  to  make  as  data  for  the  high  Reynolds  numbers  cannot  be  ob¬ 
tained  in  laboratory  facilities,  when  Froude  numbers  and  draft 
ratios  are  to  be  reproduced  correctly  as  well. 
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Postcritical  CD~values  can  be  obtained  for  smooth  circular  cylin¬ 
ders  at  subcritical  Reynolds  numbers  by  applying  trip  wires  on  ine 
cylinder  surface.  The  present  model  tests  have  revealed  that 

stable  force  coefficients  are  obtained  on  a  1:50  scale  model  of  a 
ship  equipped  with  bilge  keels  and  trip  wires  for  the  range  of 
Froude  and  Reynolds  numbers  tested.  The  CD~value  was  found  to  in¬ 
crease  for  a  1:25  scale  model.  This  increase  may  be  related  to  the 
decrease  in  turbulence  level  of  the  ambient  flow  for  the  tests  1 

with  the  1:25  scale  model.  Another  and  more  important  observation  , 

was  made  during  testing.  When  the  ship  was  free  to  roll,  pitch  and 
heave  the  drag  force  increased  by  as  much  as  40%  compared  to  the 
value  found  when  the  ship  was  restrained  at  the  still  water 
drafts,  with  zero  heel.  The  increase  is  due  to  the  heel  and  sink- 
age  of  the  ship  caused  by  the  pressures  acting  on  the  ship  hull. 

This  observation  may  be  part  of  the  reason  for  the  substantial 
scatter  in  CD~values  reported  in  the  literature.  Furthermore,  the 
effect  has  obviously  important  consequences  for  the  design  of  the 
mooring  system. 

Finally,  tests  have  been  conducted  to  investigate  if  instability  , 

motions  of  a  violent  character  would  occur  within  the  range  of  ■ 

Froude  numbers  from  0.08  to  0.16.  No  severe  motions  that  could  ' 

lead  to  capsizing  or  grounding  of  the  ship  were  observed.  However, 
the  static  heel  and  sinkage  of  the  ship  increased  with  increasing 
Froude  number.  Superimposed  on  these  average  movements  were  high  ’ 

frequency  heave,  roll  and  pitch  motions.  During  testing  it  was  ! 

noticed  that  the  roll  motion  initiated  the  pitch  and  heave  mo-  \ 

tions.  Furthermore,  it  was  observed  that  a  certain  heel  was  ne¬ 
cessary  to  have  sustained  roll  motions. 


3. 


SCOPE  OF  WORK 


The  model  study  was  planned  in  accordance  with  the  description  of 
work  for  Phase  I  of  "Investigations  and  Tests  to  Determine  Hydro- 
dynamic  Forces  and  Moments  on  Adjacent  Ships  Moored  in  a  Current". 
Phase  I  covers: 

Preliminary  tests  to  define  the  range  of  applicability  of  main 
tests . 

The  phase  I  tests  shall  be  short  yet  valid  quick-look  investi¬ 
gations  aimed  at  getting  critical  guideline  data. 

For  the  floating  body  (a  ship  model  or  an  adjusted  horizontal  cy¬ 
linder)  the  following  investigations  were  to  be  made: 

o  Find  the  maximum  flume  width  at  which  the  lateral  force  on  the 
transversely  mounted  vessel  is  not  affected  by  more  than  ap¬ 
proximately  5%  relative  to  the  situation  with  an  infinitely 
wide  flume. 

o  Determine  the  influence  of  the  Reynolds  number  on  the  lateral 
drag  force  on  the  hull  shaped  body  in  a  beam  current  so  that 
the  Phase  II  measurements  can  be  properly  converted  to  full 
scale  predictions  on  drag  forces  and  moments. 

o  Investigate  the  behaviour  of  a  floating  body  when  it  is  moored 
in  a  beam  current  and  the  speed  of  the  current  is  increased. 
Specifically  determine  if  and  when  erratic  or  violent  motions 
take  place. 

For  the  specified  types  of  investigation,  Statement  of  Work,  Ref. 
1,  specifically  describes  the  environmental  conditions  under  which 
the  tests  have  to  be  carried  out. 
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4.  MODEL  SET-UP  AND  INSTRUMENTATION 

4.1  Model  Basin 

The  set-up  was  made  in  a  rectangular  tank  with  the  dimensions  25  m 
x  35  m  where  the  inlet  and  outlet  weirs  were  located  along  the 
short  sides.  Positions  for  movable  walls  applied  in  the  blockage 
tests  were  marked  on  the  horizontal  floor  of  the  tank  as  indicated 
in  Fig.  No.  4.1  showing  the  location  of  the  walls  and  the  vessels 
relative  to  the  tank. 
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Fig.  4.1  The  Position  of  the  Vessel  and  the  Movable  Walls  in  the 
Basin. 

The  flow  through  the  tank  was  established  by  use  of  a  pump-  and 
head  tank  system  in  connection  with  inlet-outlet  flumes  around  the 
tank.  The  principle  for  the  system  operation  is  shown  in  Fig.  4.2. 


across  the  flume  in  the  level  mentioned  above.  The  horizontal  di¬ 
stribution  of  the  current  was  then  adjusted  by  lowering  or  eleva¬ 
ting  the  outlet  weirs.  When  the  uniform  horizontal  distribution 
was  established  the  vertical  velocity  profile  upstream  was  measu¬ 
red  for  control. 

In  order  to  determine  the  roughness  regime  for  the  model  flow  a 
generalized  Colebrooks  and  Whites  formula  has  been  adopted  for  the 
description  of  the  model  flow.  A  high  degree  of  correlation  has 
been  found  using  the  following  physical  parameters: 

k  =  0.001  m  Bottom  roughess. 

6  =  0.002  m  Boundary  layer  thickness. 

(Dependent  on  the  average  velocity) . 

E  =  3 0  *  1 0  J  Reynolds  Number  based  on  the  depth. 

V  =  0.13  m/s  Average  velocity  (3  ft/sec.  (prototype). 

The  calculated  as  well  as  the  measured  velocity  profile  are  shown 
in  Fig.  4.3.  The  agreement  between  the  two  is  good. 
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easured  and  Calculated  Velocity  Profile. 


>s  used  in  this  study  were  scale  models  (1:50  and  1:25) 
prototype  ship  with  the  principal  dimensions  as  out- 
)le  4.1.  The  body  lines  of  the  vessel  are  shown  in  Fig. 


e  modeJ 


Dimension  Prototype  1:50  1:25 


Lpp  (m) 

160.90 

3.218 

6.437 

T  (m) 

8.0 

0.  160 

0.320 

A  (m3 ) 

18541 

0.144 

1.154 

Bm  (m) 

23.0 

0.460 

0.920 

Bilge  radius  (m) 

3.05 

0.061 

0 . 122 

Heights  of  decks 

(m  above  base  line) 

Main  at  FP 

16.764 

0.335 

0.671 

Main  at  AP 

15.24 

0.305 

0.610 

Main  at  mid 

13.716 

0.274 

0.549 

Second  at  FP 

12.497 

0.250 

0.500 

Second  at  AP 

12.802 

0.256 

0.512 

GM  (m) 

2.134 

0.043 

— 

Ballasted  to 

0-trim 

Table  4.1  Dimensions  for  Vessels.  L  :  length  between  perpendicu¬ 
lars.  T:  draft  of  vessel.  A:  Displacement.  Bm:  moulded 
breadth.  FP:  fore  perpendicular.  AP:  aft  perpendicular. 
GM:  transverse  metacenter  height. 

During  the  tests  (excluding  some  Reynolds  number  tests)  the  model 
ship  was  equipped  with  bilge  keels.  The  height  of  the  keels  cor¬ 
responded  to  0.45  m  and  covered  the  central  66.0  m  of  the  vessel. 
The  vessels  were  tested  with  rudder  and  propeller. 

During  some  Reynolds'  number  tests  the  models  were  equipped  with 
trip-wires  in  order  to  stimulate  turbulence  in  the  boundary  layer 
flow.  These  wires  extended  from  the  ends  of  the  bilge  keels  to  the 
perpendiculars.  The  precise  locations  of  the  trip  wires  can  be 
seen  from  Fig.  No.  4.4.  The  diameter  of  these  wires  corresponds  to 
8*10  ^  times  the  bilge  radius  which  should  trig  the  turbulence 
most  efficiently  according  to  Ref.  3. 


The  moorings  were  simulated  by  stiff  rods  irrespective  that  natu¬ 
ral  moorings  are  elastic.  The  large  stiffness  of  the  mooring  rods 
is  also  responsible  for  large  fluctuations  in  the  mooring  forces 
as  the  energy  absorption  of  the  stiff  moorings  is  minimum. 

In  addition,  the  horizontal  movements,  surge,  sway  and  yaw  were 
prevented  by  the  use  of  stiff  mooring  rods. 

The  model  moorings  were  anchored  on  the  main  deck  at  the  axis  of 
the  ship  at  both  perpendiculars.  The  heights  above  the  keel  level 
are  given  in  Table  4.1. 

Anchored  barges  were  simulated  by  fixed  anchoring  points  100  m 
upstream  of  the  vessel.  The  anchoring  points  were  elevated  corres¬ 
ponding  to  1.5  m  above  the  water  level  giving  the  moorings  an  in¬ 
clination  of  approximately  2.5  degrees. 

At  the  fixed  mooring  points  and  at  the  vessel  the  moorings  were 
allowed  to  rotate  by  use  of  spherical  bearings.  Thereby  the  action 
of  moments  caused  by  the  single  moorings  was  prevented. 


Fig.  4.5  Mooring  Configuration. 


4.5  Vessel  Suspension 

In  order  to  measure  the  effects  of  pure  blocking  and  to  study 
changes  in  forces  due  to  changes  in  the  Reynolds  number  the  ves¬ 
sels  were  suspended  from  thin  wires  during  these  tests  keeping  the 
bottom  of  the  vessels  horizontal.  The  drag  forces  would  therefore 
not  be  influenced  by  reduction  of  the  underkeel  clearance  other¬ 
wise  appearing  during  heave,  roll  and  pitch  movements. 

On  the  other  hand  the  differences  in  the  flow  around  the  bow  and 
the  stern  may  introduce  vertical  movements,  and  these  may  become 
larger  in  a  narrow  flume.  The  blocking  effects  can  therefore  not 
be  separated  from  natural  effects  in  an  infinitely  wide  flume. 
Thus,  the  suspension  chosen  is  a  compromise  in  order  to  investi¬ 
gate  the  effects  due  to  pure  blocking  of  the  flume. 


During  the  instability  tests  the  vessel  was  floating  and  thus  al¬ 
lowed  to  have  the  vertical  motions,  whereas  the  horizontal  motions 
were  prevented  by  the  rigid  mooring  arrangement. 


4.6  Measuring  Techniques 

Current 

The  current  was  measured  by  means  of  either  a  bi-directional  ultra¬ 
sonic  current  meter  giving  the  instantaneous  speed  continuously  or 
by  a  propeller  current  meter.  The  output  from  the  former  device  is 
an  analogue  electrical  voltage  proportional  to  the  velocity  and 
from  the  latter  a  number  on  a  counter  preset  to  a  certain  measur¬ 
ing  period. 

During  the  instability  tests  and  the  Reynolds'  Number  tests,  the 
current  speed  was  measured  indirectly  by  means  of  a  circular  plate 
mounted  on  a  force  transducer.  The  drag  coefficient  on  a  circular 
plate  is  independant  of  the  Reynolds  Number.  The  theoretical  value 
of  the  drag  coefficient  is  1.19,  a  value  that  has  been  verified 
experimentally  in  a  number  of  studies. 

The  current  force  acting  on  the  circular  plate  has  been  used  as  a 
reference  during  the  instability  tests.  The  plate  was  mounted  up¬ 
stream  of  the  vessels  in  a  depth  corresponding  to  the  level  of  the 
forces  acting  on  the  vessels. 

M°2Ei22_I!2E222 

The  mooring  rods  were  fastened  to  force  transducers  consisting  of 
a  beam  equipped  with  a  pair  of  strain  gauges.  The  output  is  a  vol¬ 
tage  proportional  to  the  deformation  of  the  beam  and  thereby  pro¬ 
portional  to  the  acting  force. 


Water  Level 


Resistance  type  wave  gauges  consisting  of  two  parallel  rods  were 
used  to  determine  the  position  of  the  water  level  in  the  flume 
during  testing.  In  addition  to  this,  a  mechanical  point  gauge  was 
used  as  reference. 


Movements 

The  vertical  movements  of  the  (floating)  vessel  were  measured  by 
angular  potentiometers  at  the  bow  and  mid  ship  starboard  and  port 
side. 

The  pure  motions  heave,  roll  and  pitch  are  easily  derived  by  use 
of  summation/ subtraction  amplifiers. 


2§!:§_Analysis 

The  electrical  signals  from  the  various  gauges  are  fed  into  a  32- 
channel  microcomputer  for  online  analysis.  This  analysis  gave  as 
result  for  each  channel: 

the  absolute  maximum 
the  absolute  minimum 
the  average  value 
the  root  mean  square  value 

for  the  total  time  series  and  for  five  subseries. 

Selected  time  series  were  continuously  recorded  on  a  12-channel 
UV-recorder.  From  these  records  it  was  possible  to  obtain  an  over¬ 
all  impression  of  the  behaviour  of  the  forces  acting  on  the  moored 
vessel  and  its  motions. 

During  the  instability  tests  the  time  series  of  the  motions  and  of 
the  forces  were  recorded  on  a  7-channel  analogue  tape  recorder. 
These  signals  were  later  digitized  and  stored  on  a  microcomputer 
making  it  possible  to  produce  plots  of  the  time  series  of  the  pure 
motions  and  of  the  resultant  forces. 


Estimated  Inaccuracies 


Table  4.2  below  outlines  the  various  sources  of  inaccuracies  and 
their  estimated  contributions  of  the  uncertainties  of  the  total 
force  (or  force  coefficient) . 


Item 

Absolute 

Contribution 
to  calculated  C„ 

- E7 - 

Vessel 

Length 

1  mm 

-4 

3*10 

Draft  (ballasting) 

1 . 5  mm 

9  *  2  0 “3 

Water 

Depth 

1.5  mm 

6- 10'3 

Average  current  speed 

Level  of  turbulence 

5  mm/  s 

8* 10~2 

( uncon trolable) 

5- 10-2 

Eguigment 

Current  meter 

~  5  mm/s 

8 • 10~2 

Force  Transducer 

5  grammes/kg 

5- 10-3 

Total  relative  standard 

deviation  on 

-2 

force  coefficients 

8-10 

Notation 


The  lateral  drag  force  coefficient  for  the  models  has  been  deter¬ 
mined  as: 


CD  *  V<*  P  LPP 


V>) 


where  is  the  lateral  drag  force  coefficient,  F^  the  total  late¬ 


ral  force,  i.e.  the  sum  of  the  forces  measured  fore  and  aft. 


is 


the  water  density,  L  the  length  between  perpendiculars  and  T  the 

P  P 

draft.  V  is  the  velocity  measured  at  a  depth  corresponding  to  0.53 


times  the  draft.  This  velocity  deviates  insignificantly  from  the 
velocity  averaged  over  the  draft  of  the  ship.  The  current  meter 
was  mounted  at  this  depth  during  the  tests. 
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BLOCKING  TESTS 


5.1  Problem  Assessment 


It  is  a  well  known  fact  that  the  forces  experienced  by  structures 
placed  in  a  confined  flow  are  larger  than  those  in  an  infinite 
flow  area.  The  phenomenon  has  mainly  been  studied  for  2-dimensio- 
nal  flow  around  cylindrical  structures,  for  which  analytical  or 
semi-empirical  correction  formulas  have  been  proposed.  The  drag  on 
a  cylinder  in  confined  flow  may  be  found  from  the  expression: 


c  =  c 

D  DM 


(1  + 


+  7  E  CDM>'2 


(5.1) 


where  e  is  the  blockage  ratio,  CDM  the  measured  drag  coefficient 
and  CD  the  corrected  value  of  this. 

For  a  floating  ship  in  shallow  waters  exposed  to  a  beam  current 
the  flow  field  is  much  more  complicated  and  no  correction  formulas 
are  readily  available.  Applying  the  above  Equation  (5.1)  indicates 
that  for  a  CDM  value  of  1.5  and  a  correction  being  less  than  5%, 
the  value  of  e  should  be  less  than  1/15,  i.e.  for  a  draft  ratio  of 
0.6  the  width  of  the  flume  should  be  9  times  the  length  of  the 
ship.  However,  a  major  part  of  the  flow  passes  underneath  the  ves¬ 
sel,  and  hence  a  much  narrower  flume  is  expected  to  give  values 
within  5%  difference  from  the  infinite  area  case. 


5.2  Model  Tests 

The  blocking  tests  were  carried  out  using  the  1:50  model  ship.  The 
range  of  flume  width  to  the  ship  length  ratio  was  varied  from  1.5 
through  2.0,  3.0,  3.5,  4.0  to  7.8  with  the  Froude  number  0.08  in 
order  to  establish  a  reference  case.  For  the  Froude  numbers  of 
0.12  and  0.16  tests  were  made  with  two  blockage  ratios  of  3.0  and 


4.0.  All  tests  were  conducted  with  a  draft  to  depth  ratio  of  0.6. 


It  was  observed  that  a  floating  ship  heeled  several  degrees  at  the 
Frouae  number  of  0.08,  for  which  reason  the  vessel  was  kept  fixed 
in  the  vertial  direction  in  order  to  prevent  the  roll,  heave  and 
pitch  motions  as  described  in  Section  4.  It  was  also  observed  du¬ 
ring  the  testing  that  the  floating  and  heeling  ship  experienced 
larger  transverse  forces  than  the  fixed  ship.  Further  discussion 
of  this  is  given  in  Section  6. 


5.3  Test  Results 

The  resulting  drag  coefficients  for  the  situation  with  ff  =  0.08 
are  given  in  the  table  below.  The  values  are  average  values  over 
25  min.  test  time  corresponding  to  approximately  3  hours  in  na¬ 
ture  . 


B/L 

V  (m/s) 

(Prototype) 

CD 

1.5 

0.788 

2.44 

2.0 

0.846 

1.71 

3.0 

0.754 

1.61 

3.5 

0.780 

1.51 

3.5 

0.803 

1.48 

4.0 

0.745 

1.64 

7.8 

0.791 

1.32 

Table  5.1  Drag  Coefficients,  I  =  0.08.  Vertically  Fixed  Ship. 


Tests  were  made  with  Froude  numbers  0.12  and  0.16  at  the  two  flume 
widths  B/L  ~  3.0  and  4.0.  Table  5.2  below  summarizes  the  results 
with  respect  to  drag  coefficients. 
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>#• : 
V«V 

(»ViI 


VJ 


V  (m/s) 


1.46 

1.28 

1.37 

1.36 


1.71 

1.73 

1.76 

1.75 


1  .70 
1.84 

1.57 

1.58 


1.67 

1.65 

1.49 

1.55 


Table  5.2  Drag  Coefficients,  I  ~  0.12  and  0.16. 

The  effect  of  the  flume  width  on  the  lateral  drag  coefficient  is 
quite  clear  from  the  figures  given  in  the  tables  above. 


5.4  Discussion  of  Results 

The  drag  coefficients  found  in  the  model  tests  are  shown  in  Fig. 
5.1  below.  Apart  from  the  C^-value  for  E  =  0.08  and  B/L  =  4  the 
variation  of  CD  with  the  B/L-value  is  rather  consistent.  Taking 
the  CD~value  found  for  B/L  =  7.8  as  applicable  for  the  infinite 
wide  flume  case  it  appears  that  a  flume  width  of  4  to  5  times  the 
length  of  the  ship  will  yield  lateral  drag  coefficients  that  dif¬ 
fer  less  than  5  per  cent  from  the  infinite  flume  case. 
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LATERAL  DRAG  COEFFICIENT. 
FIXED  MODELS 


\\ 

w 

w 

\ 

V  \ 

\  \ 

\  \ 

\  V 

x 

4  N  o-* 


A  F  -  o  08 
x  F  ~  0  1  2 
•  F  ~  0  1  6 


flume  wi  dth  b/l 


Fig.  5.1  CD~Values  vs.  Blockage  Ratio,  B/L,  for  3  F-Values. 

The  lateral  drag  coefficient  as  recommended  by  OCIMF,  Ref.  /!/,  is 
shown  versus  the  depth  to  draft  ratio  in  Fig.  5.2.  As  indicated  in 
the  figure  a  somewhat  larger  CD~value  is  given  by  OCIMF  for  H/D 
=  1.67  than  found  in  the  present  study.  No  information  is  given  in 
Ref.  / 1  /  with  respect  to  blockage  ratio  in  the  model  tests  or 
Froude  number,  but  the  Reynolds  number  was  of  the  order  of  10 5  and 
it  seems  likely  that  the  ship  has  been  held  fixed  in  vertical  po¬ 
sition  as  for  the  present  investigations. 
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1  ""2 


7  H/T 


Fig.  5.2  Lateral  Drag  Coefficients,  CD,  on  a  VLCC  as  Recommended 
by  OCIMF  for  Various  Depth/Draft  Ratios,  H/T. 

The  discrepancy  between  the  OCIMF  and  the  present  data  is,  how¬ 
ever,  not  surprising  considering  that  different  types  of  ships 
were  modelled  and  different  facilities  applied.  The  reported  data 
in  literature  show  a  large  scatter  in  lateral  drag  coefficients, 
as  illustrated  by  the  data  compiled  by  NCEL ,  Ref.  / 2/,  and  repro¬ 
duced  in  Fig.  5.3. 


Relative  Drag  Coefficients  vs.  the  Depth/Draft  Ratio. 
From  Ref.  / 2 / . 


6. 


REYNOLDS  NUMBER  EFFECTS 


6.1  Assessment  of  the  Problem 

As  for  the  blockage  phenomenon,  the  effect  of  Reynols  number  on 
fluid  drag  forces  on  bluff  bodies  has  mainly  been  studied  for  2- 
dimensional  flow  around  cylindrical  structures.  The  well-known 
variation  of  CD  for  a  circular  cylinder  is  shown  in  Fig.  6.1  below 
together  with  sketches  of  the  boundary  layer  and  separation  cha¬ 
racteristics.  As  seen  from  Fig.  6.1  the  value  of  CD  is  closely 
connected  to  the  characteristics  of  the  boundary  layer  and  the 
position  of  the  separation  point.  Other  phenomena  that  affect  the 
Cp-value  and  its  variation  with  the  Reynolds  number  are  surface 
roughness  and  turbulence  in  the  oncoming  flow.  For  cylinders  in 
close  proximity  to  a  rigid  boundary,  the  drag  coefficient  is  af¬ 
fected  by  the  magnitude  of  the  gap  between  the  cylinder  and  the 
boundary.  Furthermore,  the  boundary  layer  developing  along  the 
fixed  boundary  influences  the  forces  on  the  cylinder. 

Extensive  model  test  programmes  have  been  conducted  to  investigate 
the  possibilities  of  reproducing  high  Reynolds  number  forces  on 
cylinders  at  lower  Re-values  by  some  artificial  means.  The  use  of 
roughness  or  trip  wires  to  stimulate  the  development  of  turbulent 
boundary  layers  is  the  most  frequently  applied  methods.  Ref.  / 3 / 
summarizes  the  results  of  a  recent  investigation  on  this  phenome¬ 
non,  and  it  is  concluded  that  the  use  of  trip  wires  and  Re  ~  10 5 
yielded  forces  and  vortex  shedding  frequencies  comparative  to 
those  for  smooth  cylinders  at  high  Reynolds  number  (~ 1 0 7 )  .  The 
results  are  reproduced  in  Fig.  6.2,  where  the  Cp-values  and 
Strouhal  numbers,  S,  are  given. 


Reynolds  Number 


CHANGES  FROM  SU3CRIT1CAL  (1.2.31  TO  POSTCRIOCAL 
181  FLO V/;  NATURAL  SEQUENCE  OF  CHANGES  FOR  A 

smooth  cylinder  (diagrammatic; 

Fig.  6.1  CD  vs.  the  Reynolds  Number  together  with  Sketches  of 
the  Flow  and  the  Boundary  Layer.  Ref.  / 3 / . 
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I  5  90- tO 
o0r  :4,o° 


SIMULATED  CONDITIONS-' 
TRIP  WIRES v 


•NATURAL  CONDITIONS 


3  90-10 
:  50°  v 


0  0  5m  CYLINDER  I  0  25m  CYLINOER  | 


0  2m  CYLINOER 
'  IN  No  3  W  T 


la)  DRAG  COEFFICIENT 


lb)  VORTEX- SHEDOING  TRIP  WIRES 

FREQUENCY  3  9D-102 


NATURAL  CONDITIONS  - 


,0*  2  2  4  5  6  8  }0S  2  3^56  •  ltf6  2  3  ^50  ®  10  J 

COMPOSITE  plot  TO  SHOW  THE  DECREE  OF  INVARIANCE  ACHIEVED  WITH  TRIP  WIRES 

Fig.  6.2  Effect  of  Trip  Wires  on  Drag  Coefficient,  CD,  and 
Strouhal  Number,  S.  From  Ref.  / 3 / . 

For  bodies  with  more  sharp-edged  corners  the  drag  coefficient  be¬ 
comes  more  constant  and  for  plates  CD  is  in  fact  invariant  for 
Reynolds  numbers  above  10 3 . 

A  ship  hull  fitted  with  bilge  keels  has  features  which  resembles 
both  the  circular  cylinder  and  a  more  sharp-edged  body,  meaning 
that  the  Reynolds  number  effects  a  priori  are  expected  to  be  less 
significant  than  for  the  circular  cylinder. 

It  is  important  to  note  that  for  flow  around  ships  the  Froude  num¬ 
ber  may  have  significant  influence  on  the  force.  Therefore  the 
Reynolds  number  effects  cannot  be  investigated  just  by  increasing 
the  velocity  as  this  also  increases  the  Froude  number.  Only  by 
using  two  different  models  at  different  scales  the  influence  of 


% 
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Reynolds  numbers  can  be  investigted  for  otherwise  identical  condi¬ 
tions  . 


6.2  Model  Tests 

The  existence  of  possible  Reynolds  number  effects  on  the  lateral 
current  force  acting  on  model  ships  has  been  investigated  through 
a  test  series.  Two  different  scale  models  of  the  same  prototype 
ship  have  been  used  for  this  purpose.  The  two  models  (of  scale 
1:25  and  1:50)  have  been  exposed  to  beam  current  of  approximately 
identical  Froude  numbers  but  different  Reynolds  numbers.  In  both 
situations  the  models  were  in  vertically  fixed  positions  and  the 
lateral  drag  coefficients  have  been  determined  and  compared. 

In  order  to  obtain  an  accurate  comparison  between  the  forces  on 
the  two  ships,  the  force  excerted  on  a  circular  plate  was  used  as 
reference.  The  plate  was  placed  at  a  certain  relative  depth  below 
the  surface  and  distance  upstream  the  ships. 

The  key  figures  from  the  tests  are  then  the  ratio  between  the  cur¬ 
rent  force  per  unit  area  on  the  vessel  and  on  the  circular  plate. 

The  flume  width  for  these  tests  was  chosen  as  small  as  2  times  Lpp 
making  it  possible  to  obtain  a  reasonable  Froude  number  for  the 
large  scale  model  (1:25).  Although  this  B/Lpp  ratio  is  known  to 
yield  significant  blockage  effects,  it  has  no  importance  for  the 
relative  force  ratios  for  the  two  ships. 

The  large  model  was  tested  in  two  conditions:  With  and  without 
bilge  keels  on  the  hull.  The  small  model  has  also  been  investiga¬ 
ted  with  and  without  the  bilge  keels,  and  furthermore  tests  with 
both  bilge  keels  and  trip  wire  were  made. 

The  trip  wire  is  applied  to  stimulate  a  turbulent  boundary  layer 
on  the  hull  sections  outside  the  bilge  keels. 


6.3  Test  Results 

For  the  small  model  the  following  results  were  obtained: 


Froude 

Number 

V/v'gD 

Reynolds 

Number 

V* B  /v 
m 

Turbulence 

Stimulators 

Relative 

Drag  Coeffi¬ 
cient 

F  /F 

Y/  C 

0.16 

1.04*  105 

Bare  hull 

1.10 

0.12 

7.84* 104 

Bare  hull 

1.05 

0.08 

5.21- 104 

Bare  hull 

1.19 

0.16 

1.04-  105 

Bilge  keels 

0.94 

0.12 

7.84-  104 

Bilge  keels 

0.93 

0.08 

5.21*  104 

Bilge  keels 

1.05 

0.16 

1.04- 105 

Bilge+trip  wire 

0.94 

0,12 

7.84-104 

Bilge+trip  wire 

1.00 

0.08 

5.21-  104 

Bilge+trip  wire 

0.95 

Table  6.1  Force  Ratios  for  Model  Vessel.  Scale  1:50. 


It  is  seen  that  the  situation  with  bilge  keels  is  very  similar  to 
the  situation  with  both  bilge  keels  and  trip  wires,  and  further 
that  the  situation  with  bare  hull  deviates  significantly  from  the 
two  other  situations.  The  drag  coefficients  for  the  bare  hull  are 
greater  than  for  the  hull  equipped  with  turbulence  stimulation 
devices . 

For  the  large  (1:25)  model,  the  ratio  F^/F^  was  found  to  equal 
1.29  for  naked  hull  as  well  as  for  the  case  with  bilge  keels  fit¬ 
ted.  The  Froude  number  was  approximately  0.06  and  the  Reynolds 
number  1 . 1  •  10  5 . 

The  test  results  are  depicted  in  Fig.  6.3  as  function  of  Froude 
and  Reynolds  number. 


Fy  /  Fc 


o  bare  hull 
&  B  I  LGE  KEELS 
0  BILGE  KEELS ' 
TRIP  WIRES 


SCALE 

model 


4  I  SCALE  MODEL 
F  SO  08 


Re  •  10 -5 


(0.08) 


Fig.  6.3  Relative  Force  Coefficients  found  in  the  Model  Tests, 


6.4  Discussion  of  Results 

The  influence  of  the  Reynolds  number  on  the  lateral  force  coeffi¬ 
cient  has  not  been  reported  extensively  in  the  literature.  One 
reason  for  this  is  naturally  the  difficulties  inherent  in  achiev¬ 
ing  prototype  values  of  the  Reynolds  number  without  otherwise  di¬ 
sturbing  or  changing  the  flow  picture.  In  water  tanks  limitations 
to  flow  speed  as  well  as  model  size  are  present.  One  should  also 
bear  in  mind  that  especially  for  shallow  water  problems  Frouae 
number  effects  interfere  with  the  phenomenon  when  large  current 
speeds  are  applied.  Wind  tunnel  testing  puts  even  more  severe  re¬ 
strictions  on  model  size,  but  larger  flow  velocity  can  be  achieved 
and  tests  applicable  for  water  flows  at  relatively  small  Frouae 
numbers  can  be  conducted.  It  is,  however,  quite  evident  that  full 
scale  Reynolds  numbers  cannot  be  achieved  in  any  of  the  two  types 
of  facilities.  The  only  feasible  avenue  for  obtaining  exact  values 
of  prototype  forces  on  vessels  is  to  conduct  full  scale  experi¬ 
ments.  The  problems  associated  with  such  experiments  are  evident, 
and  full  scale  data  are  very  sparse. 
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In  Ref.  / 4/  Treshchevsky  and  Korotkin  have  reported  results  of  a 
prototype  experiments  and  comparisons  with  wind  tunnel  model 
tests.  Based  on  information  about  the  influence  of  roughness  on 
the  drag  coefficient  of  circular  cylinders,  roughness  of  a  certain 
height  was  glued  to  the  hull  of  the  model  ships  and  the  resulting 
lateral  force  coefficients  compared  remarkably  well  with  the  model 
prototype  data  as  shown  in  Fig.  6.4  below.  The  smooth  hull  model 
yielded  somewhat  smaller  coefficient  values.  This  may  be  due  to  a 
transition  in  the  boundary  layer  which  occurs  on  cylinders  at  Rey¬ 
nolds  number  of  2*  10  s.  The  details  about  the  experiments  are  not 
given  in  Ref.  / 4 / ,  e.g.  the  depth  to  draft  ratio  and  the  value  of 
the  Froude  number  are  not  given. 
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Fig.  6.4  Lateral  Force  Coefficient  on  Full  Scale  Ship  and  Two 
Models.  From  Ref.  / 4 / . 

English  and  Wise,  Ref.  /  5  /  ,  have  also  investigated  the  influence 
of  turbulence  stimulating  devices  on  the  lateral  drag  force  on 
model  ships.  Instead  of  distributed  roughness  they  applied  trip 
wires  on  the  model,  a  method  which  a  recent  investigation  (see 


29 


dhi 


Ref.  /3/)  on  cylinder  drag  has  shown  to  be  most  appropriate  in 
simulating  high  Reynolds  number  flow.  A  decrease  of  7%  in  CQ  was 
found  when  trip  wires  were  fitted  to  the  naked  hull.  When  bilge 
keels  were  added  CD  increased  by  27%.  This  increase  was  followed 
by  changes  in  the  flow  picture  around  the  hull  caused  by  the  se¬ 
paration  zone  induced  at  the  bilge  keels.  Again  some  details  about 
the  test  parameters  are  missing,  e.g.  Reynolds  number  value  and 
the  depth  to  draft  ratio.  The  author  arrives  at  the  conclusion 
that  scale  effects  {Re  effects)  are  minimized  by  applying  bilge 
keels  on  the  model  supplemented  with  trip  wires  extending  from  the 
extremities  of  the  bilge  keels  to  the  bow  and  stern. 

As  shown  previously  in  Fig.  5.2  the  OCIMF,  Ref.  / 1 / ,  has  reported 
Cp-values  for  various  depth/draft  ratios  for  VLCC ' s  based  on  model 
tests.  Again,  details  of  the  tests  are  not  given  in  Ref.  /!/,  but 
it  is  indicated  that  the  Reynolds  numbers  in  the  tests  varied 
between  1.2  and  1.8  times  10  s,  and  that  a  decrease  in  Cp-values 
was  noted  with  increasing  Re  number  within  this  interval,  and  that 
changes  in  the  flow  around  the  bow  and  stern  were  associated  with 
the  decrease  in  CD~value:  Vortex  shedding  disappeared  at  the  lar¬ 
ger  Re-numbers.  The  form  of  the  VLCC '  s  has  similarities  with 
(smooth)  circular  cylinders,  and  the  decrease  in  CD  may  have  rea¬ 
sons  as  for  the  cylinder  case.  The  magnitude  of  the  decrease  is 
not  given,  but  the  resulting  CD~values  should  be  those  associated 
with  the  largest  Re-numbers.  By  analogy  to  the  circular  cylinder 
this  may  seem  somewhat  non-conservative,  as  vortices  are  re-estab¬ 
lished  at  large  Re-numbers,  where  CQ  increases  above  the  value 
found  at  critical  Re-numbers  of  2  to  5  times  10 5 . 

The  NCEL  report  by  Palo  and  Owens,  Ref.  /2 /,  points  to  an  even 
larger  uncertainty  with  respect  to  the  CD~values  and  the  influence 
of  Reynold  number  on  these.  Fig.  6.5  below  from  Ref.  / 2 /  shows  the 
C^-values  found  from  various  calculation  models  and  recommenda¬ 
tions.  The  scatter  in  the  Cp-values  is  significant. 


Ktfurc  Nondi mcnatonii  Ittcrtl  force  coelf icarncv 


Fig.  6.5  Lateral  Drag  Forces  on  Ships  in  Steady  Current.  From 

Ref.  / 2 / . 

The  data  obtained  in  the  present  project  indicate  that  the  appli¬ 
cation  of  bilge  keels  and  trip  wires  on  the  hull  of  the  1:50  scale 
vessels  leads  to  stable  force  coefficients  for  the  range  of  Rey¬ 
nolds  and  Froude  numbers  tested.  This  is  in  agreement  of  findings 
presented  in  Refs.  /  4  /  and  /  5  / .  The  reasons  for  the  increase  in 
the  lateral  force  ratio  for  the  1:25  vessel  are  not  obvious.  For 
the  naked  hull  case  the  increase  is  8%  and  for  bilge  keels  fitted 
it  is  23%.  Compared  to  the  smaller  vessel  the  Reynolds  number  is 
more  than  two  times  larger,  going  from  0.5  to  1.1  times  10  s 
(Froude  number  for  the  1:25  vessel  was  somewhat  below  0.08,  appro¬ 
ximately  0.06).  There  is  no  reason  for  an  increase  in  drag  coeffi¬ 
cient  due  to  the  change  in  Reynolds  number.  During  testing  it  was 
observed  that  the  flow  was  more  stable  for  the  1:25  scale  model 
tests,  with  less  turbulence  of  the  oncoming  flow.  For  circular 
cylinders  turbulence  is  known  to  reduce  the  drag  coefficient,  when 
the  Reynolds  number  is  below  the  critical  range  beginning  at  2- 
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10 5 ,  hence  this  effect  may  partly  explain  the  difference  found  in 
the  present  tests.  Although  it  is  evident  that  models  with  bilge 
keels  (and  trip  wires)  fitted  show  more  stable  force  coefficients 
than  naked  hull  models  it  may  be  too  optimistic  to  conclude  that 
prototype  values  of  CD  are  obtained  with  these  models  at  much 
lower  Re  numbers.  However ,  a  more  sound  basis  for  drawing  conclu¬ 
sions  will  require  either  full  scale  investigations  or  a  more  ex¬ 
tensive  model  test  programme,  presumably  conducted  in  a  large  wind 
tunnel.  Such  tests  were  beyond  the  scope  for  the  preliminary  and 
limited  investigations  reported  here.  A  further  discussion  of  more 
extensive  model  tests  is  included  in  the  second  part  of  this  re¬ 
port  "Phase  II  Model  Test  Programme". 

The  results  obtained  indicate  that  turbulence  may  have  a  signifi¬ 
cant  effect  on  the  lateral  forces.  An  even  more  pronounced  effect 
was  obtained  by  changing  the  set-up  from  the  one  with  the  model 
suspended  in  wires  to  the  one  with  the  ship  free  to  have  vertical 
motions  applied  for  the  instability  tests.  The  forces  were  increa¬ 
sed  by  as  much  as  40%  in  the  latter  set-up.  This  indicates  that 
due  considerations  must  be  given  to  this  effect  when  designing  the 
prototype  mooring  arrangements. 


7. 


INSTABILITY  TESTS 


7.1  Assessment  of  the  Problem 

Oscillations  of  flexible  or  elastically  suspended  structures  im¬ 
mersed  in  a  moving  fluid  may  either  be  caused  by  external  fluctua¬ 
ting  forces,  e.g.  due  to  turbulence  in  the  fluid,  or  it  may  be  a 
result  of  the  interaction  between  the  fluid  and  the  structure.  For 
this  latter  case  two  types  of  fluctuating  forces  may  appear.  One 
type  which  is  present  even  when  the  structure  is  kept  in  fixed 
position,  e.g.  vortex  shedding  from  bluff  bodies,  and  another  type 
which  is  generated  due  to  an  initial  motion  of  the  body.  This  lat¬ 
ter  type  of  oscillating  forces  and  motions  is  traditionally  refer¬ 
red  to  as  being  of  the  instability  type. 

Large  scale  turbulence  in  natural  currents  may  contain  consider¬ 
able  energies  at  frequencies  close  to  the  natural  frequency  of  the 
vessel  motions.  An  investigation  of  this  phenomenon  has  not  been 
included  in  the  present  project. 

For  a  ship  in  shallow  water  exposed  to  a  beam  current  the  two  ty¬ 
pes  of  forces  caused  by  the  fluid-structure  interaction  may  be 
present  and  a  brief  outline  of  these  are  given  below. 

The  flow  around  the  ship  is  sketched  in  Fig.  7.1.  In  the  regions 
at  the  bow  and  stern  vortices  are  being  shed  from  the  hull  without 
flow  re-attachment  as  depicted  in  Fig.  7.1.  a  and  b.  Time  varying 
pressure  fields  acting  on  the  hull  at  the  bow  and  stern  are  re¬ 
lated  to  these  pressure  fluctuations.  If  the  resulting  forces  are 
out  of  phase  a  pitch  motion  may  be  initiated. 
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Fig.  7.1  Sketch  of  the  Flow  around  a  Ship  exposed  to  a  Beam 
Current  in  Shallow  Waters. 

In  the  central  part  of  the  hull  where  bilge  keels  are  fitted  the 
flow  field  is  sketched  in  Fig.  7.1. c.  At  the  upstream  water  sur¬ 
face  a  separation  zone  develops  with  small  return  flow  velocities. 
At  the  upstream  bilge  keels  flow  separation  occurs,  followed  by 
re-attachment  and  finally  flow  separation  occurs  at  the  downstream 
bilge  keel.  Behind  the  ship  a  large  roller  is  created  with  a  typi¬ 
cal  horizontal  extension  equal  to  the  beam  of  the  ship,  meaning 
that  the  velocity  at  the  water  surface  is  in  the  upstream  direc¬ 
tion,  as  indicated  in  Fig.  7. I. a  and  c. 

The  pressure  field  acting  on  the  hull  results  in  a  moment  which 
tends  to  give  the  ship  a  mean  roll  position.  If  the  ship  heel  is 
increased  (see  the  sketch  in  Fig.  7.2)  the  flow  area  is  decreased 
at  the  upstream  bilge  and  increased  at  the  downstream  keel.  As  the 
flow  discharge  underneath  the  ship  cannot  change  instantaneously 
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Fig.  7.2  Sketch  of  the  Flow  of  the  Bilge  Keel  Section  and  Effect 
of  Roll  Motion  on  Velocity  and  Pressure  Field. 

this  will  result  in  increased  and  decreased  flow  velocities  at  the 
up-  and  downstream  keel,  respectively,  and  associated  with  this  is 
changes  in  the  pressures  leading  to  an  amplification  of  the  turn¬ 
ing  moment  in  the  direction  of  the  roll  motion:  The  mechanism 
creating  a  traditional  instability  type  of  motion  is  thus  present. 


7.2  Model  Tests 

For  these  tests  the  ship  was  free  to  have  pitch,  roll  and  heave 
motions  whereas  the  yaw,  surge  and  sway  motions  were  restrained  by 
the  rigid  rods  applied  for  the  mooring.  During  the  tests  the  ship 
motions  and  forces  were  recorded  together  with  the  force  on  a  cir¬ 
cular  plate  yielding  a  reference  force.  The  measuring  system  is 

shown  in  Fig .  7.3. 


latory  roll  motion.  The  double-amplitude  of  this  motion  is  not 
significantly  influenced  by  the  Froude  number. 


The  heave  motion  has  a  similar  trend  with  the  Froude  number.  The 
mean  value  increased  with  increaseing  ff-values,  the  amplitudes  of 
the  oscillating  part  are  not  influenced  by  the  Froude  number. 

Finally,  the  pitch  motions  also  appear  to  have  oscillation  ampli¬ 
tudes  independent  of  the  Froude  number. 

The  motions  of  the  ship  have  been  documented  by  a  video  recording. 


7.4  Discussion  of  Results 

The  tests  have  revealed  that  time-varying  motions  of  the  moored 
ship  are  likely  to  occur  for  all  flow  conditions.  The  part  of  the 
motions,  which  may  be  due  to  vortex  shedding  is  also  expected  to 
appear  at  full  scale  as  vortex  shedding  is  recognized  to  be  pre¬ 
sent  also  for  very  large  Reynolds  numbers. 

Although  it  was  not  attempted  to  prevent  one  of  the  possible  mo¬ 
tions  to  see  the  effect  on  the  others,  the  general  impression  from 
the  observations  made  during  the  testing  was,  that  the  roll  motion 
was  the  one  initiating  the  pitch  and  the  oscillating  heave  mo¬ 
tions.  Furthermore,  it  also  seemed  that  some  heeling  was  necessary 
to  initiate  the  time-varying  roll  motion. 


Finally,  it  is  important  to  note  that  the  drag  force  increased 
substantially  when  the  vertical  restraints  on  the  ship  were  re¬ 
moved.  An  increase  of  40%  in  the  CD~value  was  found  for  one  case. 


Besides  the  apparent  consequence  for  the  design  of  the  mooring 
system  this  observation  also  influences  the  model  test  programme, 
which  is  outlined  in  the  second  part  of  the  report. 
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